arXiv:hep-lat/0503021v2 18 Jul 2006 


hep-lat/0503021 


Casimir Scaling in 
SU{2) Lattice Gauge Theory 


Carlo Piccioni 


Department of Physics 
New York University 
4 Washington Place 
New York NY 10003, USA 


Abstract 

By employing the multilevel algorithm in numerical Monte Carlo simu¬ 
lations, we evaluate the static potential in four dimensional SU{2) lattice 
gauge theory with no dynamical fermions, for static sources in the j = 1) | 

representations. We hnd data supporting the Casimir scaling hypothesis. 
With the same technique we show that the ratio of the logarithm of Wilson 
loops in the j = 2 and j = ^ representations, as a function of the area of the 
loops, also satisfies Casimir scaling. 
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1 Introduction 


Numerical evaluations of Wilson loops are of key importance in the lattice 
approach to the conhnement problem, because of the relation between Wilson 
loops and the energy between static sources. Wilson loops become increas¬ 
ingly small as the separation of the sources increases. As a result they are 
hard to measure, because statistical fluctuations cover their signal. When 
studying conhnement this is clearly a problem, because we are interested in 
the large distance properties of the sources. For this reason it is common 
to hnd examples in the literature where new algorithms are introduced, in 
order to either improve the signal, see for example the smearing [1], varia¬ 
tional [2,3] and blocking [4] algorithms, or reduce the statistical huctuations, 
see for example the multihit [5] and multilevel [6] algorithms. In particular 
the multilevel algorithm has been shown to be a very powerful error reduction 
technique to evaluate small magnitude Wilson loops [6-9]. When consider¬ 
ing Casimir scaling in lattice gauge theory [10,11], the signal is much weaker 
because higher dimensional representation Wilson loops should be obtained 
by exponentiation of fundamental loops with the ratio of the corresponding 
Casimir operators. The Casimir scaling hypothesis precisely states that the 
ratio between the static potential in some representation of the gauge group 
and the static potential in the fundamental representation at the same dis¬ 
tance, is equal to the ratio of the corresponding Casimir operators. This 
should be true from the onset of confinement to the onset of color screen¬ 
ing. Using the language of the /c-string studies performed in [12-14], at 
asymptotically large distances string tensions should be dependent on the 
Wality k of the representation of the sources rather than the representation 
itself, and as consequence a string is expected to decay into the correspond¬ 
ing stable A;-string. There is considerable evidence of Casimir scaling in the 
non-perturbative regime from Monte Carlo calculations. Stndies of the SU (3) 
case can be found for example in [15,16], while the four and three dimensional 
SU{2) cases are discussed in [2,10] and [11,17,18] respectively. Nevertheless 
there is no clear theoretical explanation, in fact Casimir scaling is known to 
be true in lattice and continunm four dimensional Yang-Mills theories only 
perturbatively. One may explain Casimir scaling in four dimensions by di¬ 
mensional rednction arguments [10], remembering that it is true exactly in 
two dimensions. However, as pointed ont in [19], in two dimensions there 
is no color screening and therefore the asymptotic behavior of Wilson loops 
at large distances is different. Other models have been proposed, based on 
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arguments on vortex theory [20] and diffusion of Wilson loops [19]. Evidence 
of exponentiation of Wilson loops with Casimir operators can also be found 
in [21], 

In this work we apply the multilevel algorithm to the study of Casimir 
scaling in four dimensional SU{2) lattice gauge theory. In this case, as in 
SU{3), the only non trivial string tension at large distances is the fundamen¬ 
tal, because we can only have k = 0,1 corresponding to the trivial and the 
fundamental representations respectively. Nevertheless it is still interesting 
to study the static potential of sources in higher dimensional representations 
at intermediate distances, where Casimir scaling is conjectured. For SU{2) 
the ratio of the Casimir operators increases very rapidly with increasing di¬ 
mension of the representation, and the determination of static potentials is 
rather difficult. We test the Casimir scaling hypothesis for the j = 1, |, 2 

representations. In section 2 we introduce static potentials and Casimir scal¬ 
ing, in section 3 we explain the technique and show the results, in section 4 
we conclude with some remarks. 


2 Casimir scaling of static potentials 


The static potential is the energy Eq of the ground state of a pair of heavy 
static sources. In the continuum formulation one can show that its relation 
to the Euclidean expectation value of Wilson loops is given by 

OO 

W{R,T) = (T(/2)|e-^'^|T(i?)) = ^ (1) 

n=0 


where |4/(i?)) is the state describing the pair at distance R and H is the 
Hamiltonian of the system with orthonormal eigenstates |<h"'(i?)), and eigen¬ 
values En{R)- The state |T(i?)) is not an eigenstate of the Hamiltonian. We 
can extract information on the ground state energy Eo{R) = V{R) by taking 
the limit T ^ oo where the sum is dominated by the first term, obtaining 

W{R,T) = Co{R)e-^^^'>^, ( 2 ) 


where Co{R) = |(T(i?)|<h°(i?))p is the overlap of |T(i?)) with the ground 
state |<h°(i?)). From this equation we obtain the static potential 


V{R) 


= — lim 

T—5-00 


\nW{R, T) 
T 


( 3 ) 
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On a lattice of spacing a Eqs.(2,3) hold with R = Ra, T = Ta. 

We work with the SU{2) lattice Yang-Mills theory in the quenched ap¬ 
proximation, described by the lattice action 

S(U)=fiY,{l-RrU^) (4) 

n 


where Un is the fundamental plaquette, P = ^ and go is the lattice coupling 

9o 

constant. We consider the Euclidean expectation values W{R,T) of Wilson 
loops WciU) = ^Tr{Y[i(zc^i) around rectangular loops C of spatial exten¬ 
sion R and temporal extension T. It is possible to consider sources in the 
higher dimensional representations Dp j = 0, 1, |, 2,of SU{2) and the 

corresponding Wilson loops obtained replacing each link of the loop with its 
representation. From the dehnition of representation and the character of 
SU{2) representations we have 


WiiU) 


2j + l 


-Tr 


D^{U.) 


1 sin[(2j + l)f 
2j + 1 sin(f) 


( 5 ) 


where Ui = Yliec^i ^ obtained expanding the element Ui G SU{2) 
in terms of the Pauli matrices as Ui = cos(|)J 2 + f sin(|)h ■ a. This formula 
allows us to calculate Wp{U) once we hnd the parameter p of Ui. From 

Eqs. (2,3), written in terms of Wj{R,T) = {Wp{U)), we can evaluate the 
static potential Vj{R) between static sources in the representation j. 

According to the Casimir scaling hypothesis, [10,11], the ratio between 
the static potential in some representation and the static potential in the 
fundamental representation at the same distance, is equal to the ratio of the 
corresponding Casimir operators 


V,(R) Cl 

2 2 


( 6 ) 


where Cj = j{j + 1) is the Casimir operator of the representation j. If the 
static potential is linearly rising with the distance and the Casimir scaling 
hypothesis holds true, we have Vj{R) = ajR, where aj is the string tension in 
the representation j, and therefore aj = |j(i + From Eq. (3) we have 

for Wilson loops in the representation j, In Wj {R, T) = -|-1) In Wi (R, T) 
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at large T, from which we hnd 


Wj(R,T) 


T) 




( 7 ) 


We see immediately that the Wilson loops decrease dramatically as j in¬ 
creases. For example if a Wilson loop in the fundamental representation 
is of the order of 10“^, then in the j = 1, |,2 representations should be of 
the order of 10“®, 10“^°, 10“^® respectively. To measure Wilson loops of these 
small magnitudes, it is necessary to introduce algorithms that either enhance 
their signal or decrease the statistical fluctuations. This is the subject of the 
following section. 


3 Technique and results 


On the lattice we can write Eq.(2) for the representation j in the form 
Wj{R,T) = , where we have introduced the static potential 

in lattice units Vj{R) = Vj{Ra)a. Then the static potential in lattice units 
is given by 


V,(R) 


— lim In 

T—>oo 


Wj{R,T + l) 
W,{R,f) 


( 8 ) 


In practice for hxed R we valuate the quantities Vj{R, T) = — In 

and extrapolate Vj{R) when a plateau value in T is reached. We use lattices 
of size 16^ and 12^ x 16 for j = \ and j = 1 respectively, and lattices of size 
8^ X 12 for i = I and j = 2, at P = 2.5. We update the gauge helds by use of 
the heat bath algorithm introduced in [22]. To enhance the projection on the 
ground state we use spatial smeared links [1] in the Wilson loops, obtained 
by the following prescription 


Wj{R,T+l) 

Wj{R,f) 


where 


U^{n) 


Ui{n) + eSj 
^J<let{Ui{n) + eSi)' 


(9) 


= Y,[U,{n)UPn + (n + ?) + U] (n - - j + ^)], (10) 

and = 1,2,3. The procedure is iterated a number n* of times. We hnd 
that the choices of the smearing parameters e = 0.1, ns = 37 for j = 
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Figure 1: (color online). The physical static potential for the j = 1, | 

representations. The dots are the numerical values obtained, and the solid 
lines are the potential hts. 


e = 0.3, = 13 for j = I and e = 0.5, = 31 for j = \ work well. 

For j = 2 it is difficult to determine the smearing parameters, in fact the 
signal is very weak and the dependence of the parameters on the spatial 
extent of the loops is more sensitive. On the temporal links of the loops, 
we apply the multilevel algorithm as an error reduction technique. This 
algorithm was introduced in [6] as generalization of the multihit algorithm [5], 
where each link is substituted by its expectation value, obtained keeping the 
neighboring links hxed. We tried to apply the multihit algorithm with no 
satisfactory results. The multilevel algorithm is employed with a two level 
scheme according to which, following the notation and terminology of [6], the 
Euclidean expectation values W{R,T) = {Wc{U)) are given by 

W{k,f) = (L(0)„,{l|T(0)]lT(l)|]...llT(f-2)]|T(t-l)|]}„^,jL(r)J(), (11) 

where L denotes the spatial transporters, T denotes the two-link operators 
and [.] denotes the time-slices expectation values, obtained keeping hxed the 
links on the boundary of the time-slices. With this scheme only even temporal 
extension Wilson loops can be calculated. We evaluate the [T] averages from 
7 updates of the time-slices of width 1, and the [[T][T]] averages from 9 
updates of the time-slices of width 2. Note that since the Wilson loops are 
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Figure 2: (color online). Left: ratio of the static potentials for j = 1 and 
j = |, against the ratio of the corresponding Casimir operators. Right: ratio 

of the logarithm of Wilson loops for j = 2, against ^ = 8. 
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obtained with the expectation values of the two-link operators T, we cannot 
make use of Eq. (5) to evaluate the higher dimensional representation traces, 
and we must consider the matrix form of the links in the representation j. 
In figure 1 we plot the results obtained for the physical static potentials for 
j = 1, |. We obtain the following least square fits for the potentials in 

lattice units 

h(«) = b - I + (12) 

it 

with the following values for the parameters 

= 0.54(1) ei = 0.24(1) <ji = 0.035(1) 

Vfl = 1.40(3) Cl = 0.62(3) <ti = 0.094(5) 

= 2.4(4) e| = 1.0(3) as = 0.15(8). 

The lattice size a = 0.087(l)/m, or a = 0.440(6)GeI4“^, is extrapolated 

0-1 

from the relation cri = -f-, giving to the string tension the physical value 
^ /oT = 425MeI4. In figure 2 we plot the ratios of the static potentials 

V 2 Vi{H) 

in lattice units, against the ratios of the Casimir operators ^ = | for j = 
1, and 7 ^ = 5 for j = |. For j = 2 we plot the ratios as a 

’Cl -2 2 -J ^ InlVi(A) 

function of the area, against ^ = 8. In this case evaluation of the static 

potential is difficult since the signal is very weak and mostly lost in statistical 
fluctuations. Exponentiating the fundamental potential rescaled with ratio of 
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the Casimir operators, we find that Wilson loops in the j = 2 representation 
should get as small as ~ 10“^^. 

4 Conclusions and outlook 

The data obtained support the Casimir scaling hypothesis in four dimen¬ 
sional SU (2) lattice gauge theory. By employing the multilevel algorithm we 
are able to extract the static potential from Wilson loops, for the j = 1, | 

representations. In these cases we test directly the Casimir scaling hypoth¬ 
esis, by considering the ratio of the lattice potentials with the fundamental. 
In the j = 2 representation case we obtain an indirect test of Casimir scal¬ 
ing by considering the ratio of the logarithms of Wilson loops as a function 
of the area. Analysis of the data highlights the multilevel algorithm as a 
powerful technique to evaluate small magnitude Wilson loops. Although the 
string picture would be more appropriate to a larger length scale than the 
one studied in this work, below 0.7/m, we could not consider larger lattices 
due to the limited computational availability. 

As observed in the introduction, Casimir scaling is expected to be true 
up to a certain distance between the sources after which the gluonic string 
holding the pair together should break or be screened to the fundamental, 
for integer or half integer j respectively. Adjoint string breaking in four 
dimensional SU{2) lattice gauge theory has been observed in [2,23]. In 
[23] the mass of the resulting gluelump is measured and the string breaking 
scale is set to Rb = 1.25/m. A rough estimate of the decay scale for the 
j = I representation is set at | of the adjoint in [14]. Since our interquarks 
separations are much below these decay scales, we cannot test color screening 
and the A^-ality dependent string tension scenario. This would certainly be 
possible using larger size lattices, resulting in new interesting numerical data 
that, although trivially expected on theoretical grounds, are not yet available 
in the literature. 
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